Heat conduction in simple networks: Controlling heat flow through inter-chain coupling 
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The heat conduction in simple networks consisting of different one dimensional nonlinear chains is studied. 
We find that the coupling between chains has different function in heat conduction compared with that in electric 
current. This might find application in controlling heat flow in complex networks. 
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Energy and information transports on networks, such as the 
metabolism network, neuronal network, porous material net- 
work, and oil production network, etc., have been studied for 
a long time [1] and is recently getting more attention because 
of the hectic activity in the complex networks and the great 
progress in nanoscale fabrication technology where the nao- 
tube/nanowire networks can be made for different purposes 
|8y. It is found that the electric transport 
changes linearly with the number of added bonds IHQl. The 
whole resistance of network can be figured out by the Kirch- 
hoff second law for the complicated parallel and serial electric 
circuit fl. 

However, little is known about heat conduction in the com- 
plex networks, although some progress has been achieved in 
the study of heat conduction in single one dimensional chains 
(See Ref.|@] and the references therein). The fundamental 
question for heat conduction in one dimensional chains is that 
what is the necessary and/or sufficient condition for the heat 
conduction to obey the Fourier law. From computer sim- 
ulations, it is found that in ID nonlinear lattices with on- 
site potential such as the Frenkel-Kontorova (FK) model and 
the <fi 4 model, the heat conduction obeys the Fourier's law, 
namely, the heat conductivity is size independent] 10], which 
is also called normal heat conduction. Whereas in other non- 
linear lattices without on-site potential, thus momentum is 
conserved, such as the Fermi-Pasta-Ulam (FPU) and alike 
models, the heat conduction exhibits anomalous behavior jl ill , 
namely the heat conductivity k diverges with the system size 
N as k ~ N s . A great effort has been devoted to under- 
stand the physical origin and the value of the divergent ex- 
ponent <5 lll2ll . It is found that the anomalous heat conduction 
is due to the anomalous diffusion and a quantitative connec- 
tion between them has been established] _13fl . Most recently, 
we found that both the normal and anomalous heat conduc- 
tion can be described by an effective phonon theory under the 
same framework] 14]. 

More importantly, it is found that a single ID chain con- 
sists of two different lattices exhibit very interesting physical 
phenomena such as thermal rectification! 15] and negative dif- 
ferential thermal resistance llflL Experimental on nanotube 
has verified the rectification 11711 . Opening up a new field of 
controlling heat flow from simple (or complex) nano scale net- 



works. 

All studies on the ID single chain can be regarded as the 
first step to understand the heat conduction on realistic situa- 
tions, i.e., complex networks. In general, a complex network 
consists of many ID (or quasi- ID) chains with a diversity of 
couplings among them. Therefore, the key to understand the 
heat conduction on networks is to understand the influence 
of coupling to heat fluxes in simple networks, i.e., coupled 
chains. According to the best of our knowledge, this problem 
has not been investigated so far. 

For the sake of simplicity, we would like to consider m 
ID chains with several couplings between any two of them. 
To be more specific, we take the FPU-/? chain Hill as the 
basic element and each chain are contacted with the Nose- 
Hoover thermostat Hill at the two ends, keeping the first 
and the last particle of the chain at temperature 7), and Tj, 
respectively. Without coupling, each chain has a Hamilto- 
nian H = Y^ihPi + v ( x h x i+i)> where V(xi,x i+1 ) = 
\ Sifci+i - x i) + I J2i( x i+i ~ x i) 4 > x i represents the dis- 
placement from the equilibrium position of the i'th particle. 
The motion of the particles for i = 2, 3, ■ • ■ , N— 1 satisfy the 
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equations for the heat baths are = ^ — 1, = — 1. 
The dynamical equations for the first and last particles are 

The temperature is defined as T(i) = (pf) and the heat 
flux along the chain is J = (Pt gf^ )- Suppose there is a 
coupling between the node i of one chain and the node j of 
another chain, then we have an additional new potential V-j = 
\{Xi — Xj) 2 + \{xi — Xj) 4 . The equations of node i and node 
j become 
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To investigate the influence of coupling by numerical sim- 
ulations, we take Th = 0.7 and 7] = 0.5 for all the chains 
and first consider the case of m = 2, i.e., two coupled chains 
in this Letter. The two chains are coupled at different nodes 
We find that both the temperature distribution and the 
total flux in the steady state are changed with the coupling 
positions. 

Case I: two chains of the same length coupled together 
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FIG. 1: Temperature distributions of two coupled chains of length 
N — 20 with different coupling positions. The thin lines (green) 
denote the coupling. The insets are the schematic configurations of 
coupled chains, and the arrows there label the direction of heat flow, 
(a) No coupling, (b) Coupling at i = j = 10. (c) Coupling at i = 
10, j = 15. (d) Coupling ati = 5, j = 15. 
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FIG. 2: The corresponding fluxes of Fig. Q] 



The two chains are identical of length TV = 20. Like 
all models of heat conduction, there are always temperature 
jumps at the two boundaries lfl9ll as is clearly shown in Fig. [TJ 
When two chains are coupled together regardless of the cou- 
pling position, there is also temperature jump at the junction, 
seeFig.QIb)-(d). 

In Fig. |2] where we shows the corresponding fluxes of Fig. 
Q]where the arrows denote the directions of fluxes. Fig. [2ja) is 
easy to understand from their identity, where two uncoupled 
chains have the same flux. 

Fig. |2fb) shows a very interesting result - the reduction 
of the heat current. This is completely different from elec- 
tric circuit. It is well known that a circuit of two chains with 



four equal resistance R connected by a conduction line at the 
middle is a symmetric circuit. Since there is no potential dif- 
ference between the two connecting points, there is no current 
through the middle connection line, thus the current in the cir- 
cuit does not change! It remains the same if the two chains 
are disconnected. 

What makes the "thermal circuit" different from the elec- 
tric circuit? To this end, we need to go to the definition of 
temperature. The temperature is a measure of the kinetics 
of the particle. It is an ensemble (time) average of the ki- 
netic energy. Without coupling, the middle particle at each 
chain is connected only by its two nearest neighbors. After 
coupling, the middle particle is connected with three particles 
which changes its equation of motion. Even though the two 
particles in the middle have the same temperature (same aver- 
age kinetic energy and same velocity distribution), it does not 
mean that the two particles always oscillate in the same way. 
This is the fundamental difference between the electric circuit 
and thermal circuit. 

In fact, the coupling of the second chain to the first chain is 
equivalent to the introduction of an interface resistance at the 
junction. This resistance is also called the Kapitza resistance 
which is defined as: Rint = AT/ J, where the AT is the 
temperature jump between the left and right particles of the 
interface (coupled particle in the middle). Therefore the heat 
current through each chain is: 



J 



T h -Ti 
2R + Ri n 
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which is obviously less than Jo = Tfe 2 ^ T ' for the uncoupled 
chain. 

In the case of without any coupling, the temperature of the 
i'th particle inside the FPU chain is: 



T ~ T h -\AT h \ 



i - 1 
N-2 



(T.-IAT.I-T-IATQ, (3) 



where |AT,| and | AT; | is the temperature jump at the both 
ends between the heat bath and the first/last particle of the 
chain, respectively. The heat current flows at the junction can 
be understood from this formula. For instance, the particle at 
i = 10, has higher temperature than the particle of i = 15. 
Heat flows always from high temperature to low temperature, 
therefore, if one connects i = 10 in upper chain to particle 
i = 15 in lower chain, there will be heat current flows from 
i = 10 (higher temperature) in upper chain to particle i = 
15 (low temperature) in the lower chain. This will drag the 
temperature of particle i = 10 down a little bit, thus we see 
the increase of the heat current in the part of i E [2, 10] in 
upper chain in Fig. 2(c) compared with the case in Fig. 2(b). 
In contrast, as the heat current flows to particle i = 15 at 
lower chain, the temperature at i = 15 is increased, thus the 
increase of the temperature difference between i = 15 and 
i = 20, which leads to the increase of heat current in segment 
of i G [15, 20] in lower chain. This is what we observe in Fig. 
[T](c). The same mechanism applies also to Fig[TJ(d). 
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Case II: Two chains of different length coupled together 

We would like to extend above ideas to more general case, 
namely, coupling of two chains of different length and multi- 
ple coupling. We find that the reduction of heat flux by cou- 
pling is quite general. Fig. [3] shows the temperature distribu- 
tion of two coupled chains with different lengthes Ni = 20 
and N2 = 30, respectively. (See Figure caption for more in- 
formation.) It is easy to see that Fig. [3] has some similarity 
with Fig. Q] i.e., there are temperature jumps at the coupled 
particles and the coupled particles have the approximate same 
temperature. 

Another interesting thing is that the crossing couplings 
make the middle part of the coupled chains appear a tempera- 
ture plateau which might be useful in heat control. 

Fig. |4]shows the corresponding fluxes of Fig. [3] The longer 
chain, N2 = 30 has smaller heat current. Although the cou- 
pling of the two chains of different length is not at the sym- 
metrical point, there is still no current through the coupling. 

In fact, the like in the previous case shown in Fig. (1) and 
(2), the heat current flow in the (multi) coupled chain of dif- 
ferent length can be also understood from Eq. j3J. According 
to this formula, we can roughly estimate that the T5 at short 
chain is roughly the same as T§ in the longer chain. Thus there 
is no heat flow between them. The only influence is the intro- 
duction of an interface resistance which drags down the heat 
current through each chain as is seen in Fig. [3jb). 

It is not difficult to estimate that T5 at chain N = 30 is 
larger than T5 at chain N = 30. This is why we see the 
current flow from particle 5 at lower (longer) chain to particle 
5 at upper (shorter) chain in Fig. |4] However, the temperature 
of particle 15 at upper chain (shorter) is almost the same as 
the temperature of particle 25 at lower chain (longer) if the 
two chains are uncoupled. However, due to first coupling, the 
temperature of particle 15 at upper chain is slightly increased, 
this is why we see the current flows from particle 15 in upper 
chain to the particle 25 in low chain as shown in FigUc). 

More complicated and more interesting case is shown in 
Fig. Ilfd), where we have two crossing couplings: i\ = 5 
is connected to 12 = 15, and j\ = 20 is connected to 22 = 
10. Use Eq.®, we can again estimate that Tj 1= 5 > Ti 2= i$, 
thus we see the current flow from upper chain to the lower 
chain. Similarly, there is heat current flows from lower chain 
(particle j'2 = 10) to upper chain (ji = 20). 

We have also checked the heat conduction in multiple cou- 
pled chains with a diversity of couplings, such as in the three 
coupled chains with different lengthes, and observed the sim- 
ilar results as in the case of two coupled chains. We conclude 
that, in general, the coupling will introduce an interface re- 
sistance at the junction thus affect the heat flow through the 
whole system. 

Case III: Single chain with loop 

Another interesting question is how the self-coupling or a 
shortcut in a single chain affect the heat current? This case 
happens very frequently in the polymer chain and biological 
systems???. For example, if there is a shortcut between the 
node i and the node j of a chain (see the inset of Fig. Ufa)), 



0.7 T 



0.6- 



0.7 



0,6 



0.5 1 


0.7r^ 



(a) 



:o 



- 0.5 

30 
-1 0.7 



\ 

Y* \ 


/ ' T 


°°o 0o **** 

(b) 


** 

**« 4 

**** 



30 



0.6- 



0.5 



(C) 



20 



0.6 



- 0.5 
30 



* 



***** 



(d) 



30 



FIG. 3: Temperature distributions of two coupled chains with dif- 
ferent lengthes Ni = 20 and N2 = 30 with different coupling 
positions. The thin (green and purple) lines denote the coupling. 
The insets are the schematic configurations of coupling chains, and 
the arrows indicate the direction of heat flow, (a) no coupling, (b) 
one coupling added at i — 5, j = 8; (c) two couplings added at 
ii — 5, ji = 5 and 12 = 15, j'2 = 25; (d) two crossing couplings 
added at ii = 5,ji = 20 and 12 = 15, j'2 = 10. 
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FIG. 4: The corresponding fluxes of Fig. [3] 



does this shortcut reduce the flux of the chain? In the case of 
traffic flow (reference ????), the shortcut increases the capac- 
ity of traffic because the vehicles have more free space to go. 
However, in the thermal circuit, the flux should be reduced 
because of the interface resistance. The line with "stars" in 
Fig. |5jb) shows the result with N — 20, i — 5, and j = 15. 
Comparing it with Fig. [2ja) of no coupling, it is easy to see 
that the flux is reduced almost 50%. 

We also study the dependence of the reduction of the heat 
flux on the coupling strength k. From the definition of junc- 
tion resistance we know that the larger the degree of destroy- 
ing the correlation between the coupled particle and its neigh- 
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bors is, the larger Ri n t, i.e., Rint should monotonously in- 
crease with k. When the correlation is completely destroyed, 
Rint cannot be increased more by further increasing k. There- 
fore, there is a saturation effect for R int and the effect of junc- 
tion resistance when k is large enough. Let's confirm this pre- 
diction by numerical simulations. In this situation, the cou- 
pling potential becomes 

V/ j = ^(x i -x j )* + ^x i -x J )'. (4) 

Substituting Eq. (0]i into Eq. (Q]i we get the dynamical equa- 
tions for the particles with coupling strength k. Our numerical 
simulations show that for a single chain with the self-coupling, 
the larger the coupling is, the more reduction of flux. Fig. 
|3b) shows three typical cases where the lines with "circles", 
"stars" and "squares" denote the cases of k — 0.5 and 2.0, 
respectively. From the middle parts of this figure it is ease 
to see that the larger coupling makes less flux go through the 
original path. We notice that the coupling also changes the 
temperature distribution. The strong the coupling is, the two 
particles connected by the coupling have more close temper- 
atures, as shown in Fig. [2a). For observing the influence 
of coupling strength in more detail, Fig. |5jc) shows how the 
fluxes change with the coupling strength k where the line with 
"circles" denotes the total flux and the line with "stars" the 
flux going through the shortcut. Obviously, the total flux be- 
comes stabilized when k > 1 and the flux through the shortcut 
is monotonously increase with k, confirming the saturation ef- 
fect. The saturation effect has been also observed in the cou- 
pling of two coupled chains, see Fig. [2d) for how the total 
flux of the two chains in Fig. |2jb) changes with the coupling 
strength k. 

In conclusions, we have studied the influence coupling in 
simple networks on the heat conduction. It is found that dif- 
ferent from the electric circuit, the coupling affect very much 
the heat current flow in the thermal circuit. Any introduction 
of coupling is equivalent to an introduction of an interface 
resistance, thus influence largely the heat current in the cir- 
cuit. The study may shed lights for studying heat conduction 
in complex networks. 
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